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An  experimental  study  is  performed  on  a  three-bladed  rotor  model  in  two  water  tanks.  The  blade 
pitch,  rotational  velocity,  descent  angle,  and  descent  speed  are  all  varied  in  order  to  simulate  a 
wide  range  of  rotorcraft  operating  states,  focusing  on  descent  cases  where  the  rotor  is  operating 
in  or  near  vortex  ring  state  —  an  area  in  which  there  is  currently  very  little  available  data.  Flow 
visualization  is  done  by  injecting  air  bubbles  and  fluorescent  dye  tangentially  from  the  blade  tips 
to  mark  the  vortex  core,  showing  the  development  of  both  short-wave  (“sinuous”)  and  long-wave 
(“leapfrogging”)  instabilities  on  the  helical  vortices  in  the  wake.  Strain  gages  are  used  to  record 
transient  loads,  allowing  a  correlation  between  the  rotor  thrust  performance  and  the  development 
of  the  vortex  wake.  Reynolds  numbers  are  of  order  10r>  and  test  runs  are  performed  for  extended 
periods  —  up  to  500  rotor  revolutions  —  demonstrating  the  repeatability  of  the  patterns  of  thrust 
variation.  The  data  indicate  that  as  the  instabilities  develop,  adjacent  vortices  merge  and  form 
thick  vortex  rings,  particularly  during  descent.  Periodic  shedding  of  these  rings  from  the  wake 
associated  with  vortex  ring  state  is  observed,  resulting  in  peak-to-peak  thrust  fluctuations  of  up 
to  95%  of  the  mean  and  occurring  at  regular  intervals  of  20—50  rotor  revolutions,  depending  on 
flow  parameters. 


Notation 

c  blade  chord 

A  rotor  disk  area,  n R2 

Ct  thrust  coefficient,  T  /  pAVfip 

L  blade  length 

N  number  of  rotor  revolutions 

R  rotor  radius 

Rec  Reynolds  number  based  on  chord,  cVtip/v 
T  rotor  thrust 
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V  free  stream  (towing)  speed 

Vh  velocity  induced  at  rotor  in  hover,  {T/2pA)1i2 
Vt.ip  rotor  tip  speed,  flR 

Vx  rotor  forward  flight  speed,  V  cos  a 
Vz  rotor  descent  speed,  V  sin  a 

a  descent  angle 

9  collective  pitch  angle  at  0.75  R 

A  short-wave  instability  wavelength 

v  kinematic  viscosity 

p  water  density 

a  standard  deviation  of  mean  rotor  thrust 

crs  rotor  solidity  ratio 

t  vortex  ring  shedding  period 

Q  rotor  rotation  rate 
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Introduction 

An  accurate  understanding  of  the  physics  of  helical 
vortex  wakes  has  long  been  regarded  as  one  of  the  most 
difficult  problems  in  fluid  dynamics.  With  implications 
on  the  performance  of  propellers,  wind  turbines,  and 
helicopter  rotors,  the  problem  is  of  practical  interest 
to  many.  Even  before  the  days  of  modern  production 
helicopter  flight,  the  issue  of  the  nature  and  stability 
of  ring  vortices  and  helical  vortices  had  been  analyzed 
extensively.  Levy  and  Forsdyke  (Ref.  1)  performed  a 
stability  analysis  on  a  single  helical  vortex  in  1928, 
and  more  recently,  Landgrebe  (Ref.  2)  and  Widnall 
(Ref.  3)  have  added  to  and  corrected  this  study.  Gupta 
and  Loewy  (Ref.  4)  have  performed  a  similar  stability 
analysis  on  multiple  interdigitated  helical  vortices. 

Despite  the  focus  that  modern  helicopter  flight  has 
brought  to  the  problem  in  the  last  hundred  years,  and 
despite  the  power  of  modern  computers  and  experi¬ 
mental  tools,  a  true  grasp  of  the  physics  of  helical  vor¬ 
tices  has  remained  elusive.  While  somewhat  reason¬ 
able  approximations  of  their  behavior  can  be  made  un¬ 
der  highly  simplified  scenarios,  there  is  much  progress 
still  to  be  made  on  the  problem  of  real  helical  vortices. 
The  current  state  of  rotor  wake  modeling  is  that  while 
accurate  computations  of  relatively  simple  wake  flows 
(such  as  hover)  are  beginning  to  produce  results  that 
are  quantitatively  interesting,  computations  of  rotor 
wakes  that  undergo  massive  interaction  with  the  rotor 
(such  as  for  the  vortex  ring  state)  are  currently  of  only 
qualitative  interest  at  best. 

In  the  early  days  of  helicopter  flight,  a  number  of 
simplified  models  appeared  -  such  as  the  classic  mo¬ 
mentum  theory  and  the  blade  element  momentum  the¬ 
ory  (Refs.  5,6)  -  which  were  capable  of  predicting  gross 
performance  characteristics  for  a  rotor  (such  as  thrust 
and  power)  but  were  unable  to  capture  the  detailed  dy¬ 
namics  of  the  wake  flow  field.  In  view  of  the  acute  ef¬ 
fect  the  wake  has  on  the  rotor  performance,  researchers 
began  to  undertake  more  detailed  experimental  stud¬ 
ies  of  the  wake  flow  field,  beginning  with  wind  tun¬ 
nel  testing  and  smoke  flow  visualization  as  early  as 
the  1920s  (Refs.  7,  13,  14).  Later,  experiments  were 
done  using  hot  wire  anemometry  (HWA)  and  other 
techniques  (Ref.  15)  to  measure  the  velocity  field  it¬ 
self.  Even  more  recent  work  has  used  laser  Doppler 
velocimetry  (LDV)  (Refs.  16-18)  and  particle  image 
velocimetry  (PIV)  (Refs.  19,20)  as  non-invasive  means 
of  achieving  the  same. 

The  power  of  modern  computers  has  recently  been 


harnessed  by  researchers  using  more  sophisticated,  de¬ 
tailed  models  of  the  flow  field.  The  current  state  of  the 
art  in  rotor  wake  modeling  is  “free  wake  analysis”,  a 
method  in  which  the  wake  flow  field  is  calculated  based 
on  the  induced  velocities  of  all  the  tip  vortices  as  well 
as  the  inboard  vortex  sheets  (see  Refs.  21,22).  And 
while  such  techniques  have  led  to  significant  progress 
over  the  past  ten  years  in  the  understanding  of  simpler 
cases  of  flight  such  as  hover,  they  are  not  yet  capable 
of  predicting  the  behavior  of  the  wake  in  more  compli¬ 
cated  flight  states  where  the  aircraft  is  maneuvering  or 
descending.  At  this  point,  there  is  no  reliable  physical 
model  for  rotor  wake  flow,  nor  is  one  possible  without 
better  experimental  knowledge  of  the  flow  than  cur¬ 
rently  exists. 

The  problem  for  aerodynamicists  is  that  the  flight 
regimes  that  are  of  greatest  interest  to  the  rotorcraft 
community  are  precisely  the  ones  in  which  CFD  pre¬ 
dictions  are  least  accurate.  When  a  helicopter  is  de¬ 
scending  rapidly,  the  very  thing  that  makes  the  wake 
solution  so  difficult  -  the  intense  interaction  between 
the  rotor  and  its  vortex  wake  -  causes  large,  un¬ 
steady  dynamic  loads  on  the  blades.  Under  certain 
circumstances  (when  the  rotor  descent  velocity  ap¬ 
proximately  matches  the  wake  velocity),  the  aircraft 
can  encounter  a  condition  known  as  vortex  ring  state 
(VRS),  where  the  tip  vortices  merge  together,  forming 
a  thick  vortex  ring  that  remains  near  the  rotor  plane, 
disrupting  the  inflow  and  causing  a  dramatic  reduc¬ 
tion  in  lift.  This  unstable  ring  typically  undergoes 
shedding/reformation  cycles  that  result  in  large  fluc¬ 
tuations  in  thrust  that  make  the  aircraft  quite  difficult 
to  control.  As  retired  test  pilot  Mott  F.  Stanchfield 
says,  “In  my  opinion,  a  mature  VRS  is  the  most  haz¬ 
ardous  condition  that  exists  in  the  realm  of  helicopter 
aeronautics”  (Ref.  23). 

Although  the  aerodynamics  of  rotors  in  descending 
flight  -  and  in  VRS  in  particular  -  has  been  the  sub¬ 
ject  of  research  for  many  years  (see  Refs.  7-10,  24), 
that  work  has  largely  focused  on  measurements  of  the 
rotor  thrust  and  power  with  very  little  flow  field  vi¬ 
sualization,  and  even  less  simultaneous  measurement 
of  thrust /power  and  flow  fields.  This  is  likely  due 
to  the  relative  “disorder”  of  such  flows  and  the  diffi¬ 
culties  associated  with  facility/model  sizing,  turbulent 
diffusion,  and  injection  of  flow  markers.  The  present 
study  is  aimed  at  addressing  these  difficulties  and  pro¬ 
viding  a  relatively  complete  description  of  the  wake  of 
the  descending  rotor  and  its  effect  on  rotor  loading. 
Ultimately,  it  is  planned  to  develop  a  data  set  with 
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sufficient  information  to  improve  and  validate  compu¬ 
tational  methods  for  the  prediction  of  rotor  descent 
behavior. 

Flow  visualization  and  thrust  measurement  results 
are  presented  from  experiments  performed  on  a  model 
rotor  in  a  wide  range  of  operating  states  -  from  hover 
to  forward  flight  to  rapid  descent  -  with  the  focus  be¬ 
ing  on  the  vortex  ring  state  regime.  Experiments  were 
performed  in  a  60  m  long  water  tow  tank,  which  al¬ 
lowed  for  extended  run  times  that  were  more  than  suf¬ 
ficient  for  the  study  of  long-period  unsteady  flows  such 
as  those  encountered  in  descending  flight,  and  which 
were  much  longer  than  have  previously  been  performed 
in  similar  studies  (which  have  generally  been  conducted 
in  wind  tunnels  -  see  Refs.  24,25).  The  rotor’s  perfor¬ 
mance  is  quantified  by  measurements  of  its  thrust,  and 
this  information  is  correlated  with  simultaneous  flow 
visualization  images. 

The  time-history  characteristics  of  the  rotor’s  thrust 
are  examined  for  a  broad  range  of  descent  speed,  de¬ 
scent  angle,  and  collective  angle  combinations.  By  test¬ 
ing  the  rotor’s  performance  over  a  wide  variety  of  con¬ 
figurations,  the  rotor’s  performance  can  be  fully  char¬ 
acterized,  and  the  descent  conditions  in  which  VRS  be¬ 
havior  is  observed  can  be  clearly  identified.  The  thrust 
time- histories  of  these  periodic  shedding  cases  are  then 
compared  in  order  to  determine  how  the  descent  con¬ 
figuration  affects  the  amplitude,  frequency,  and  overall 
“orderliness”  of  the  observed  thrust  fluctuations.  For 
these  particular  cases,  the  flow  visualization  images  of 
the  experiment  provide  clues  as  to  the  nature  of  the 
vortex  wake  formation  and  shedding  phenomenon  that 
makes  VRS  such  a  dangerous  flight  regime. 


Experimental  Setup 

Rotor  Model 

Experiments  were  performed  with  a  three-bladed  25.4 
cm  diameter  rotor  model  featuring  manually  adjustable 
blade  pitch.  The  blades  (Fig.  1),  which  were  9.5  cm 
long,  were  molded  from  carbon  fiber  reinforced  plastic. 
The  blades  were  untapered,  with  a  1.9  cm  chord,  and 
had  a  twist  of  about  5°  (compared  with  twists  of  35- 
40°  for  typical  tilt-rotor  aircraft).  The  low  blade  twist 
was  chosen  in  order  to  minimize  rotor  separation.  The 
blade  airfoils  were  modified  ARAD-10  at  the  tip  and 
modified  ARAD-13  at  the  root.  Each  blade  had  a  0.36 


Figure  1:  Model  3-bladed  rotor.  Diameter  is  25.4 
cm.  Dye  and  air  enter  the  tubes  at  the  root,  flow 
through  the  blades  to  the  tips,  and  enter  the  flow 
tangentially. 

mm  ID  stainless  steel  tube  embedded  along  its  span, 
allowing  dye  or  air  to  be  injected  into  the  flow  from 
the  tip  to  mark  the  vortex  cores.  The  airfoil  modifi¬ 
cations  included  a  chordwise-linear  thickness  increase 
which  thickened  the  trailing  edge  and  provided  more 
room  for  the  dye/air  tubes.  In  addition,  the  aft  cam¬ 
ber  of  the  root  airfoil  was  increased. 

The  rotor  was  driven  by  a  digitally-controlled  mi¬ 
crostepper  motor  (25,000  pulses  per  revolution),  allow¬ 
ing  for  precise  control  of  both  the  rotor’s  position  and 
velocity.  The  motor  was  mounted  atop  an  89  cm  ver¬ 
tical  shaft  and  drove  a  23  cm  horizontal  shaft  onto 
which  the  rotor  was  fixed  (Fig.  2).  The  vertical  shaft 
lay  downstream  of  the  rotor  during  descent  testing,  so 
as  to  avoid  interfering  with  the  rotor’s  inflow.  Just  be¬ 
neath  the  motor  was  a  2.5  cm  thick  rectangular  mount¬ 
ing  plate  which  supported  the  model  assembly  and  also 
measured  rotor  thrust,  using  a  pair  of  flexures  instru¬ 
mented  with  120  ohm  strain  gage  bridges.  The  thrust 
readings  were  fed  to  the  computer  controlling  the  ex¬ 
periment.  Data  are  corrected  for  drag  tares  and  low- 
pass  filtered  during  post-processing  to  eliminate  high- 
frequency  electrical  and  vibrational  noise  while  retain¬ 
ing  the  main  features  of  the  signal. 

To  visualize  the  rotor’s  wake,  air  bubbles  and  sodium 
fluorescent  dye  were  released  from  the  blade  tips  in  a 
direction  tangential  to  the  blade  path.  The  dye  and 
air  were  supplied  to  the  dye  reservoir  at  the  base  of 
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the  vertical  shaft  through  thin,  flexible  plastic  tubing 
(Fig.  2).  The  dye  reservoir  was  directly  connected  to 
the  rotor  dye  tubes  through  the  horizontal  drive  shaft. 
The  pressure  deficit  at  the  blade  tips  drew  some  fluid 
from  the  supply  tube  into  the  wake,  but  in  order  to 
achieve  clear  visualization  of  the  flow  it  was  necessary 
to  supply  external  pressure. 


Figure  2:  Entire  model  assembly,  showing  locations 
of  strain  gages  and  dye  tubing.  Flow  would  be  from 
left  to  right  in  a  descent  test. 


Stationary  Tank 

Initial  testing  was  performed  in  a  1.22  x  2.44  x  1.68  m 
deep  stationary  water  tank  at  the  University  of  Cali¬ 
fornia  Berkeley.  With  the  model  fixed  in  place  atop  the 
tank,  these  tests  simulated  a  hovering  helicopter’s  flow 
field.  A  10-W  Argon  ion  laser  provided  both  planar 
and  volumetric  illumination.  For  the  two-dimensional 
lighting  tests,  a  vertical  light  sheet  was  aligned  with  the 
axis  of  the  rotor.  In  all  cases,  a  digital  video  camera 
recorded  the  flow  from  the  side  of  the  tank,  perpendic¬ 
ular  to  the  wake  direction  and  the  light  sheet.  Due  to 
the  small  size  of  the  tank,  test  runs  were  relatively  brief 
(and  between-test  intervals  relatively  long)  in  order  to 
minimize  recirculation  effects  which  would  lead  to  an 
unwanted  climb  condition. 

Generally,  air  was  used  as  the  injection  fluid  for  ini¬ 
tial  experiments  because  of  its  non-contaminating  na¬ 
ture.  Also,  in  contrast  to  the  neutrally-buoyant  dye, 
which  marked  the  vortex  cores  but  diffused  to  areas 
surrounding  the  core  as  well,  the  majority  of  the  air 
bubbles  initially  migrated  directly  to  the  low  pressure 
vortex  cores.  Thus  in  the  near-wake  of  the  rotor,  be¬ 
fore  the  vorticity  had  diffused  significantly  and  atten¬ 
uated  the  pressure  deficit  at  the  cores,  the  air  bubbles 
nicely  captured  the  details  of  the  filament  structure. 
However,  once  the  transition  to  the  far-wake  region 
occurred,  the  buoyancy  of  the  bubbles  caused  them  to 
rise  to  the  surface  quickly,  rendering  the  details  of  the 
wake  indecipherable  for  distances  greater  than  about 
one  diameter  downstream  of  the  rotor.  Later  tests  used 
neutrally-buoyant  fluorescent  dye  as  the  injection  fluid, 
which  clearly  showed  the  break-up  and  diffusion  of  the 
wake  at  greater  downstream  distances. 

In  general,  thrust  measurements  were  not  recorded 
for  the  stationary  tank  tests.  Rather,  these  tests  were 
performed  primarily  for  visualization  purposes,  as  the 
quality  of  the  images  was  significantly  better  in  the 
stationary  tank  than  in  the  towing  tank  and  the  struc¬ 
ture  and  evolution  of  the  wake  could  be  more  clearly 
discerned. 
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Steel  Frame 

(connection  to  carriage! 


Figure  3:  Side  view  of  carriage  platform  with  model 
assembly  mounted  on  turntable.  Direction  of  travel 
is  from  left  to  right. 

Towing  Tank 

The  characteristics  of  a  descending  helicopter  were  sim¬ 
ulated  by  pulling  the  model  through  water  in  the  60  m 
long  tow  tank  at  the  University  of  California’s  Rich¬ 
mond  Field  Station.  The  2.4  m  wide,  1.5  m  deep  tank 
features  a  large,  low-speed  carriage  running  along  a  set 
of  rails  on  top  of  the  tank.  The  carriage  speed,  which 
for  these  tests  ranged  from  0-50  cm/s,  could  be  con¬ 
trolled  manually  or  by  computer. 

The  model  assembly  (Fig.  2)  was  mounted  on  a  1.22 
x  1.52  m  plywood  platform,  which  was  suspended  just 
above  the  water  surface  by  a  steel  frame  connected  to 
the  carriage  (Fig.  3).  The  model  could  be  rotated 
using  a  turntable  on  the  plywood  platform,  enabling 
the  descent  angle  of  the  rotor  to  be  varied  in  0.5°  in¬ 
crements,  from  0°  (forward  flight)  to  90°  (vertical  de¬ 
scent).  A  set  of  halogen  track  lights  with  blue  dichroic 
filters  was  mounted  onto  the  front  of  the  carriage  to  il¬ 
luminate  the  flow  and  highlight  the  yellow  fluorescent 
dye.  Unfiltered  white  light  was  used  when  air  was  the 
injection  fluid.  For  visualization,  the  video  camera  was 
mounted  on  the  platform  vertically,  looking  downward 
at  the  flow  and  fixed  in  position  with  respect  to  the 
rotor. 

Results 

All  results  presented  in  this  paper,  unless  stated  oth¬ 
erwise,  refer  to  experiments  conducted  at  a  rotor  rota¬ 
tion  rate  of  fi=4  rev/s  (Vup= 319  cm/s)  and  at  a  col¬ 
lective  pitch  angle  of  0=11.6°.  This  was  taken  as  a 
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Figure  4:  Variation  in  mean  thrust  coefficient  with 
collective  pitch  angle.  Values  for  a  range  of  rotation 
rates  collapse  —  approximately  —  to  a  single  curve. 


representative  case  in  order  to  limit  the  number  of  ex¬ 
perimental  variables  and  also  because  it  produced  the 
best  flow  visualization  results.  In  addition,  since  ro¬ 
tor  thrust  is  proportional  to  the  square  of  the  rotation 
rate,  the  non-dimensional  thrust,  Ct,  should  be  inde¬ 
pendent  of  rotation  rate  (in  the  absence  of  viscous  or 
aeroelastic  effects),  being  a  function  of  the  rotor/airfoil 
geometry  and  the  blade  pitch  angle  only. 

This  result  was  tested  experimentally  for  the  case  of 
a  hovering  rotor.  The  mean  thrust  coefficient  was  mea¬ 
sured  for  a  wide  range  of  collective  angles  and  rotation 
rates  in  the  water  tow  tank.  The  expected  result  is 
that  the  individual  curves  -  marking  the  variation  in 
thrust  coefficient  with  collective  angle  -  for  each  of  the 
rotation  rates  tested  (fi=3-7  rev/s)  should  collapse  to 
a  single  curve.  However,  Fig.  4  shows  a  displacement 
of  the  five  curves  that  is  a  regular  function  of  rotation 
rate  (except  for  the  highest  speed).  For  the  rotation 
rates  tested  here,  blade  Reynolds  numbers  ranged  from 
f?ec=45, 000-106, 000,  a  region  which  is  certainly  sub¬ 
ject  to  viscous  effects.  However,  the  regularity  of  the 
displacement,  regardless  of  collective  angle,  is  surpris¬ 
ing.  The  lift  variation  seen  in  Fig.  4  could  be  accounted 
for  by  a  twist  of  about  5°,  which  is  certainly  possible 
given  the  present  blade  construction.  It  is  therefore 
considered  plausible  that  aeroelastic  twist  may  also  be 
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Figure  5:  Three-dimensional  flow  visualization  im¬ 
age  of  rotor  using  air  injection  from  tip  of  one  blade. 
Rotor  speed  is  fl=3  rev/s.  The  short-wave  sinuous 
instability  of  the  vortex  is  just  visible. 

occurring. 

The  stalling  of  the  lift  trend  at  the  highest  rotor 
speed  is  also  notable.  Whereas  the  lower  speed  curves 
show  continued  increases  in  Ct  at  collective  angles  as 
high  as  30°,  the  fi=7  rev/s  curve  plateaus  at  about  20°. 
Normally,  stall  would  be  expected  to  be  more  likely  to 
occur  at  lower  speeds.  A  possible  explanation  is  that, 
at  this  high  speed  and  collective,  the  lift  is  enough  to 
twist  the  blades  and  lead  to  a  major  separation. 

Nevertheless,  the  high  thrusts  obtained  for  all  higher 
collectives  cases  {Ct  is  over  0.03,  with  a  blade  load¬ 
ing  coefficient,  Ct/ps,  of  over  0.2)  are  surprising  and 
not  currently  understood.  Future  studies  will  include 
surface  flow  visualizations,  torque  measurements,  and 
higher  stiffness  blades  to  determine  the  nature  of  this 
behavior. 

Flow  Visualization 

Flow  visualization  testing  in  the  stationary  water  tank 
yielded  numerous  images  clearly  showing  the  develop¬ 
ment  of  the  rotor  wake  and  the  instabilities  that  cause 
it  to  break  down.  In  Fig.  5,  air  bubbles  are  injected 
from  the  tip  of  only  one  blade  for  the  sake  of  clarity. 
The  large  starting  ring  vortex  can  be  seen  on  the  left, 
expanding  and  slowing  down  as  the  rest  of  the  wake 


passes  through  it.  In  this  three-dimensional  image, 
the  short-wave  “smooth  sinuous  wave  type”  instabil¬ 
ity  -  which  was  discussed  by  Leishman  (Ref.  ?)  and 
Fukumoto  and  Miyazaki  (Ref.  11)  and  analyzed  the¬ 
oretically  by  Widnall  (Ref.  3)  and  Gupta  and  Loewy 
(Ref.  4)  -  can  clearly  be  seen  along  the  filament  in  the 
near-wake  of  the  rotor.  Ortega  et  al.  (Ref.  26)  have 
suggested  that  this  is  an  elliptic  instability  of  the  vor¬ 
tex  cores  that  develops  cooperatively  on  adjacent  helix 
turns,  though  this  observation  is  unable  to  be  verified 
here  with  air  being  released  from  only  one  blade.  In 
this  case  the  wavelength  of  the  instability  is  approxi¬ 
mately  3.75  cm,  or  2c. 

Figure  6  shows  a  series  of  two-dimensional  images 
of  the  upper  half  of  the  rotor  plane  with  dye  being 
injected  from  all  three  blade  tips.  This  view  of  the 
wake  shows  cross-sections  of  the  cores  of  the  three  tip 
vortices  and  illustrates  the  influence  that  each  vortex 
filament  has  on  its  neighbors.  The  induced  velocities 
of  adjacent  turns  cause  the  helices  to  expand  and  con¬ 
tract,  thus  altering  their  propagation  speeds  and  re¬ 
sulting  in  the  classic  “leapfrogging”  or  “vortex  pairing” 
phenomenon  often  seen  with  parallel  vortex  rings.  This 
effect  -  which  was  studied  computationally  by  Jain  and 
Conlisk  (Ref.  12)  and  experimentally  by  Ortega  et  al. 
(Ref.  26)  -  can  be  seen  in  the  pairing  of  the  second 
and  third  vortex  cores  downstream  of  the  rotor  in  (b) 
and  (c),  and  quickly  leads  to  the  complete  merger  of 
all  three  vortices  in  (e).  The  merger  of  adjacent  tip 
vortices  is  a  general  phenomenon  of  rotors  regardless 
of  the  number  of  blades  -  for  example,  Ortega  et  al. 
(Ref.  26)  observed  it  in  the  wake  of  a  two-bladed  rotor. 
However,  with  fewer  blades  the  adjacent  helix  turns 
are  farther  apart  and  thus  the  leapfrogging  and  merg¬ 
ing  processes  take  longer  to  transpire.  The  location 
where  the  three  vortices  merge  -  about  half  a  diameter 
downstream  of  the  rotor  -  can  be  taken  as  the  point 
where  the  wake  loses  its  helical  structure. 

Thrust  Measurements 

Instantaneous  thrust  measurements  were  recorded  dur¬ 
ing  the  descent  experiments  in  the  towing  tank.  These 
tests  were  conducted  over  a  range  of  towing  speeds 
from  0-50  cm/s,  descent  angles  from  a=0-90°,  and 
collective  angles  from  0=6-18°.  Descent  runs  were 
typically  performed  for  100  rotor  revolutions,  although 
some  were  conducted  for  longer  periods  in  order  to  ver¬ 
ify  the  trends  observed  during  shorter  runs.  The  data 
sampling  rate  was  200  samples  per  revolution,  and  the 
first  and  last  five  revolutions  of  the  run  were  ignored 
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Figure  6:  Two-dimensional  flow  visualization  images  of  upper  half  of  rotor  in  hover  using  fluorescent  dye 
injection  from  all  three  blade  tips.  The  “leapfrogging”  of  one  vortex  filament  over  another  can  be  seen  in 
(b)  and  (c)  as  the  three  vortices  orbit  about  each  other  and  finally  merge  in  (e). 


(due  to  the  starting  transients  from  the  carriage  mo¬ 
tion,  and  also  to  the  width  of  the  filtering  window). 

For  many  of  the  runs,  thrust  levels  remained  rel¬ 
atively  steady  over  the  duration  of  the  experiment. 
This  was  generally  the  case  for  hover,  slow  descent,  and 
very  steep  or  very  shallow  descent  angle  runs.  Figure 
7  shows  a  thrust  coefficient  time-history  for  a  hover¬ 
ing  rotor.  The  mean  thrust  coefficient  was  0.0078  and 
the  peak-to-peak  fluctuation  amplitude  was  12%  of  the 
mean.  Figure  8  shows  a  thrust  coefficient  time-history 
of  similar  form,  but  for  a  rapid  descent  at  a  fairly  shal¬ 
low  descent  angle  (a=30°).  Note  that  the  mean  thrust 
coefficient,  0.018,  is  more  than  two  times  greater  than 
in  the  hover  case,  but  the  total  fluctuation  amplitude 
is  still  only  12%  of  the  mean. 

However,  for  experiments  featuring  a  combination 
of  moderate  descent  speed  and  steep  descent  angle, 
the  rotor  thrust  characteristics  were  markedly  differ¬ 
ent.  In  some  such  cases,  the  thrust  exhibited  very  large, 
regular  fluctuations.  Figures  9-10,  for  instance,  show 
thrust  time-histories  typical  of  this  type  of  behavior. 

These  plots  exhibit  classic  VRS  characteristics,  with 
very  large,  regular  thrust  oscillations.  The  sustained 
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Figure  7:  Thrust  time-history  for  a  hover  test 
(V=0).  Mean  thrust  coefficient  is  0.0078  and  fluc¬ 
tuation  is  12%  of  the  mean. 
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figure  8:  Thrust  time-history  for  a  fast,  shallow  de¬ 
scent  run  (14/14  =  1.62,  14/14  =-0.93).  Mean  thrust 
coefficient  is  0.0181  and  fluctuation  is  12%  of  the 
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Figure  9:  Thrust  time-history  for  a  classic  VRS 
case:  a=60°,  V/14=1.25  (14/14=0.62,  14/14=- 
1.08).  Length  of  run  is  500  rotor  revolutions  instead 
of  100,  but  performance  parameters  are  unchanged. 
Mean  thrust  coefficient  is  0.0135,  fluctuation  ampli¬ 
tude  is  94%  of  the  mean,  and  fluctuation  period  is 
about  43  revolutions. 


Figure  10:  Thrust  time-history  for  another  VRS 
case  with  a=50°  and  V/14  =  1.5  (14/14=0.96, 

14/14  =-1.14).  Mean  thrust  coefficient  is  0.0173  and 
fluctuation  is  50%  of  the  mean.  Fluctuation  period 
in  this  case  is  only  about  20  revolutions. 

regularity  of  the  thrust  oscillations  was  verified  by  per¬ 
forming  long  runs  of  over  500  blade  revolutions  (see 
Fig.  9).  (For  most  runs,  a  run  of  100  revolutions  will 
suffice.) 

Other  descent  configurations  demonstrated  very  sim¬ 
ilar  oscillatory  characteristics,  but  with  significantly 
different  amplitudes  and  periods  of  fluctuation.  Fig¬ 
ure  10  shows  a  thrust  coefficient  time-history  that  ap¬ 
pears  very  similar  to  Fig.  9.  However,  in  this  case  - 
where  the  descent  angle  was  50°  and  the  towing  speed 
was  V/14= 1.5  -  the  thrust  oscillations  were  smaller 
in  magnitude  (peak-to-peak  variation  was  50%  of  the 
mean)  and  of  shorter  period  (approximately  20  revo¬ 
lutions).  The  mean  thrust  coefficient  in  this  case  was 
slightly  higher  though,  at  0.0173. 

Obviously  the  descent  configuration  greatly  influ¬ 
ences  the  oscillatory  behavior  of  the  rotor  thrust  and 
thus  can  spell  the  difference  between  a  routine  bumpy 
ride  and  a  catastrophic  loss  of  control  for  a  descending 
aircraft.  This  difference  can  be  seen  by  comparing  the 
thrust  envelopes  -  the  region  between  maximum  and 
minimum  thrust  levels  -  for  two  descent  angles  over 
a  range  of  towing  speeds.  Figure  11  shows  the  maxi¬ 
mum,  minimum,  and  mean  thrust  values  for  (a)  a=90° 


8 


0.030 


0.025 


•|  0.020  - 

0 
o 


0.015- 


\  0.010 


0.5  1.0  1.5  2.0  2.5 

Non-dimensional  Total  Speed  (V/Vh) 


0.005 


0.000 

0. 


i  r 

(b)  a  -  60°  descent 

1 

1 

< 

0 

Pre-VRS 

«—  dvrs  > 

0  0 

(aperiodic) 

Regime 
(periodic)  ^ 

-1 

+ 

0  0 

+  +  *  * 

0  0  + 

* 

»  +  +  * 

* 

Post-VRS 

0  +  +  * 

(aperiodic) 

o  0  0  0  0  0 

+  +  *  * 

I  *  ^  iti  +  ,  ,  + 
*  +  +  + 

* 

*  *  * 

*  * 

- 

* 

_ 1 _ 

0.5  1.0  1.5  2.0  2.5 

Non-dimensional  Total  Speed  (V/Vh) 


3.0 


Figure  11:  Thrust  coefficient  envelopes  —  maximum  and  minimum  thrust  levels  —  for  (a)  90°  descent  and 
(b)  60°  descent  configurations. 


and  (b)  a=60°  descent  configurations.  For  the  vertical 
descent  case  (a=90°),  the  thrust  levels  dropped  pre¬ 
cipitously  as  the  descent  speed  increased  to  about  20 
cm/s  (V/14=1.0).  However,  there  was  no  noticeable 
thrust  periodicity  accompanying  the  loss  of  lift  as  one 
would  see  with  a  classic  vortex  ring  state  case,  and  thus 
the  envelope  is  relatively  narrow.  As  the  descent  rate 
increased  further,  thrust  was  recovered  and  eventually 
reached  two  to  three  times  hover  thrust  levels  as  the 
rotor  encountered  turbulent  wake  state  and  windmill 
brake  state. 

The  rotor’s  behavior  in  a=60°  descent  is  similar, 
with  loss  of  lift  at  approximately  20  cm/s  (V/V/= 1.0), 
followed  by  recovery  and  increased  lift  as  the  descent 
speed  reaches  windmill  brake  state.  The  major  differ¬ 
ence  between  the  two  cases  is  in  the  size  of  the  enve¬ 
lope.  In  vertical  descent,  the  thrust  levels  drop  dras¬ 
tically  in  the  20-40  cm/s  descent  speed  range,  yet  the 
difference  between  maximum  and  minimum  thrust  lev¬ 
els  is  relatively  small  (compared  to  the  mean)  and  also 
somewhat  random.  However,  in  the  case  of  the  a=60° 
descent,  the  rotor  experiences  only  a  slight  reduction 
in  lift  once  VRS  is  reached,  but  also  large,  regular  os¬ 
cillations  in  thrust  as  well.  Thus  the  thrust  envelope 
for  this  case  is  noticeably  larger  -  covering  a  signifi¬ 
cant  proportion  of  the  total  thrust.  Controllability  of 
the  helicopter  under  these  conditions  is  a  major  ques¬ 


tion.  In  this  case  it  should  be  noted  that  the  variation 
in  thrust  (the  expansion  of  the  thrust  envelope)  occurs 
at  speeds  as  low  as  V/V/=  0.5  -  well  before  the  fully- 
developed  VRS  region.  The  thrust  variations  in  this 
region  are  still  aperiodic,  though  quite  large.  If  con¬ 
trolled  flight  is  sustainable  in  this  region,  this  could 
constitute  an  operationally  useful  precursor  to  the  full 
vortex  ring  state.  Note  that  no  such  precursor  is  seen 
for  vertical  descent. 

The  mean  thrust  levels  in  Fig.  11  are  also  very  re¬ 
vealing.  In  the  a=60°  VRS  region,  the  mean  thrust  is 
half-way  between  the  maximum  and  minimum  values, 
as  one  would  expect  for  a  highly-organized,  periodic 
oscillation  pattern.  But  in  the  higher  and  lower  speed 
regions  the  thrust  pattern  is  entirely  aperiodic,  and 
this  is  reflected  in  the  fact  that  the  mean  thrust  values 
are  no  longer  midway  between  the  maximum  and  min¬ 
imum  but  closer  to  the  minimum  values.  In  the  a=90° 
descent  configuration  the  mean  thrusts  appear  to  lie 
half-way  between  the  maximum  and  minimum  curves, 
but  in  this  case  the  thrust  patterns  are  aperiodic  for 
all  descent  speeds  and  the  mean  values  happen  to  lie 
in  the  middle  of  a  fairly  random  distribution  of  thrust 
values. 

Figure  12  provides  a  detailed  look  at  the  dynamics 
of  the  vortex  ring  formation  and  shedding  process  that 
precipitates  the  thrust  oscillations  discussed  above.  In 
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Figure  12:  Sequence  of  images  depicting  the  various  stages  of  vortex  ring  formation  and  shedding  during 
VRS.  Number  of  revolutions  (N)  since  the  beginning  of  the  run  is  shown  in  each  image  for  reference. 
Corresponding  thrust  coefficient  time-history  shows  the  resultant  loss  and  recovery  of  lift  during  the  cycle. 
Helicopter  at  lower  left  represents  the  rotor  orientation  and  incident  flow  direction.  Test  parameters  for 
this  case  were:  a=60°,  "V/Vh— 1.0,  12=3  rev/s,  0=11.9°  (14/14=0.50,  14/14  =-0.86).  Mean  thrust  coefficient 
was  0.0086  and  peak-to-peak  variation  was  121%  of  the  mean,  with  a  fluctuation  period  of  36  revolutions. 
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(a),  the  rotor  is  relatively  free  of  vorticity  and  the  flow 
seems  fairly  orderly  -  hence  the  thrust  is  near  its  peak 
level.  The  incoming  flow  keeps  the  tip  vortices  near 
the  rotor  plane  after  they  are  generated  and  they  be¬ 
gin  to  roll  up  into  a  thick  ring  in  (b),  interfering  with 
the  rotor’s  inflow  and  causing  the  thrust  level  to  drop 
accordingly.  In  (c),  a  strong,  tight  ring  can  be  seen 
to  have  formed  at  the  back  (downstream)  side  of  the 
rotor  as  the  thrust  approaches  its  minimum  value.  At 
this  point,  the  ring  begins  to  expand  and  separate  from 
the  back  of  the  rotor  slightly  in  (d)  and  (e),  just  as  it 
starts  to  take  form  at  the  front  (upstream)  side  of  the 
rotor.  This  “halo”  geometry,  where  the  vortex  ring  is 
approximately  equal  in  size  and  strength  at  the  front 
and  back  side  of  the  rotor  -  and  slightly  stretched  out 
as  a  result  of  the  back  side  detachment  -  roughly  co¬ 
incides  with  the  absolute  minimum  thrust  level  seen 
in  the  cycle.  Finally,  in  (f),  the  ring  has  almost  fully 
detached  from  the  rotor  plane  and  is  quickly  convected 
away  upstream  as  the  thrust  level  rebounds  sharply. 

This  entire  sequence  of  images  -  showing  the  com¬ 
plete  formation,  strengthening,  and  shedding  process 

-  spans  only  21  revolutions,  although  the  total  peak- 
to-peak  oscillation  cycle  takes  about  36  revolutions. 
Given  the  metronomic  regularity  of  the  thrust  time- 
history  for  this  run,  as  well  as  those  shown  in  Figs. 
9-10,  it  seems  clear  that  this  process  is  highly  deter¬ 
ministic  and  that  the  detachment  of  the  vortex  ring 
was  not  due  to  some  flow  disturbance  or  anomaly.  In 
fact,  the  vortex  ring  formation  process  described  above 

-  where  the  ring  takes  shape  and  gains  strength  at  the 
back  side  of  the  rotor,  then  begins  to  detach  and  stretch 
out  just  as  the  ring  starts  to  grow  at  the  front  side  - 
was  observed  during  every  thrust  oscillation  cycle  in 
this  test  case. 

Summary  Statistics 

A  primary  interest  of  the  rotorcraft  community  is  the 
determination  of  the  VRS  boundary  -  the  region  of  the 
flight  envelope  in  which  VRS  conditions  are  likely  to  be 
encountered.  Numerous  methods  have  been  employed 
for  quantifying  the  magnitude  of  the  thrust  fluctua¬ 
tion  levels  and  thus  defining  this  boundary.  Betzina 
(Ref.  24)  has  proposed  using  three  standard  deviations 
of  the  thrust  coefficient  (3<r)  normalized  by  the  mean 
value.  This  measure  is  more  consistent  than  simply 
the  peak-to-peak  amplitude  of  the  oscillations,  and  also 
less  dependent  on  the  length  of  the  run  and  will  thus 
be  used  in  this  study. 

Figure  13  shows  a  contour  plot  of  this  parameter 
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Figure  13:  Relative  magnitude  of  thrust  oscillations, 
3cr/Gra„9  (0=11.9°  runs  only). 
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Figure  14:  Amplitude  of  peak-to-peak  thrust  coeffi¬ 
cient  oscillations,  CTman;  —  Crmin  (0=11.9°  runs  only). 
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Figure  15:  Mean  thrust  coefficient  versus  non- 
dimensional  advance  speed,  14/14  (9=11.9°  runs 
only). 
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Figure  16:  Mean  thrust  coefficient  with  respect 
to  non-dimensional  advance  and  descent  speeds, 
14/17and  14/14  (0=11.9°  runs  only). 


for  all  test  runs,  plotted  with  respect  to  advance  and 
descent  speeds,  14/14,  and  14/14-  Clearly  the  area 
featuring  moderate  forward  and  descent  rates  exhibits 
the  most  dramatic  oscillations,  as  discussed  earlier.  It 
is  difficult  for  anyone  but  an  experienced  helicopter  pi¬ 
lot  to  specify  how  large  the  thrust  fluctuations  must  be 
for  the  rotor  to  be  considered  in  VRS.  However,  judg¬ 
ing  by  the  relative  magnitude  of  the  thrust  fluctuations 
shown  in  Fig.  13,  the  region  in  the  center  of  the  plot 
-  between  about  a=30°  and  70°  -  is  clearly  set  apart 
from  the  rest  of  the  flight  envelope. 

Figure  14  shows  the  peak-to-peak  amplitude  of 
thrust  fluctuations.  Although  not  apparent  from  Fig. 
13,  it  is  clear  that  the  amplitude  of  the  thrust  oscil¬ 
lations  is  still  quite  large  in  forward  flight  regimes. 
However,  the  mean  rotor  thrust  in  forward  flight  is 
significantly  greater  than  in  descent,  thus  the  relative 
magnitude  of  the  oscillations  is  less  dramatic. 

This  latter  point  can  be  illustrated  by  plotting  the 
mean  thrust  coefficient  versus  the  non-dimensional  ro¬ 
tor  advance  speed  -  the  forward  flight  speed,  14,  nor¬ 
malized  by  the  induced  velocity  at  the  rotor  plane,  14. 
This  plot,  shown  in  Fig.  15,  demonstrates  the  approx¬ 
imately  linear  increasing  trend  of  the  rotor  thrust  with 
forward  flight  speed.  This  feature  is  also  apparent  in 
Fig.  16,  a  contour  plot  of  the  mean  rotor  thrust  coef¬ 
ficient  versus  the  forward  and  descent  speeds. 

The  behavior  of  the  thrust  fluctuations  with  respect 
to  the  forward  flight  speed  is  not  quite  as  straightfor¬ 
ward,  as  shown  in  Fig.  17.  At  larger  advance  speeds 
(14/14 >1-0)  the  wake  is  quickly  swept  away  from  the 
rotor  by  the  free-stream  flow,  thus  allowing  very  little 
interaction  between  the  rotor  and  wake  and  thus  little 
fluctuation  in  the  thrust.  Likewise  at  very  low  advance 
speeds  (14/14<0.2)  the  wake  self  induction  convects 
the  rotor  away  from  the  rotor,  again  minimizing  ro¬ 
tor/wake  interaction  and  thrust  fluctuation.  However, 
at  moderate  advance  ratios  the  wake  self-induction  and 
the  freestream  flow  partially  negate  each  other,  mini¬ 
mizing  the  wake  convection  and  allowing  very  strong 
rotor/wake  interaction.  This  explains  the  wide  scatter 
of  the  fluctuations  seen  in  Fig.  17  for  14/14=0.2-1.0. 

As  discussed  above,  for  the  VRS  cases,  the  oscil¬ 
lation  frequencies  and  amplitudes  varied  significantly, 
depending  on  the  descent  configuration.  The  oscilla¬ 
tion  periods  ranged  from  20  to  50  revolutions,  and  the 
amplitudes  of  the  fluctuations  measured  from  50%  to 
95%  of  the  mean.  Figure  18  shows  the  relationship  be¬ 
tween  the  oscillation  period  and  the  non-dimensional 
advance  speed  for  the  cases  with  observable,  organized 
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Three  Sigma  /  Mean 


Figure  17:  Relative  magnitude  of  thrust  oscillations, 
■ler/C'/'  ,  versus  non-dimensional  advance  speed, 

14/14  (0=11.9°  runs  only). 


VRS-like  oscillations.  Data  are  plotted  for  a  range  of 
collective  angle  settings  (0=6-17.5°),  with  a  different 
power  law  regression  for  each  setting.  The  trend  shown 
here  is  not  unusual  -  one  would  expect  the  vortex  ring 
to  be  less  stable  and  thus  to  be  shed  more  rapidly  as 
the  advance  rate  increases.  What  is  unusual  is  how  well 
the  power  law  relationships  represent  the  data.  The 
four  regression  equations,  along  with  their  correlation 
coefficients,  are: 


0=6.0°: 

0=9.0°: 

0=14.7° 

0=17.5° 


r  =  26.6  (14/V4)-1-34 
r  =  18.6  (Vx/Vh )~1A1 
r  =  11.7  (Vx/Vft)"1-35 
r  =  9.3  (Vx/Vh)-1-70 


R2  =  0.96 
R2  =  0.98 
R2  =  0.87 
R2  =  0.79 


The  relationship  observed  here  between  oscillation 
period  and  collective  angle  -  that  for  a  given  advance 
speed,  the  shedding  period  increases  with  collective  an¬ 
gle  -  provides  some  clue  as  to  the  dynamics  of  the  vor¬ 
tex  ring  shedding  process.  As  the  collective  angle  (and, 
therefore,  the  total  lift)  are  increased,  the  strength  of 
the  tip  vortices  increases  proportionally.  In  the  vortex 
ring  state  regime,  this  equates  to  a  stronger  vortex  ring 
with  a  greater  self-induced  velocity.  This  self-induction 
is  effectively  the  force  holding  the  growing  ring  to  the 
rotor  plane  and  preventing  it  from  convecting  away 
with  the  free  stream.  Hence,  with  greater  self-induced 
(downward)  velocity  comes  longer  attachment  periods 
between  the  ring  and  the  rotor  plane. 


0  0.25  0.5  0.75  1  1.25  1.5 

Non-dimensional  Advance  Speed,  VxA/h 

Figure  18:  Variation  of  the  oscillation  period,  r, 
with  non-dimensional  advance  speed,  14/14)  for 
cases  with  observable  VRS-like  fluctuations.  The 
data  points  shown  are  the  “organized”  fluctuations 
only.  Collective  angles  range  from  0=6—17.5°. 


Conclusions 

Flow  visualization  and  thrust  measurement  exper¬ 
iments  have  been  performed  on  a  three-bladed  rotor 
model  in  two  water  tanks.  Descent  angle  and  speed 
have  been  varied  in  order  to  simulate  a  wide  range 
of  descent  configurations,  with  particular  emphasis 
on  the  vortex  ring  state  regime.  Some  of  the  specific 
conclusions  from  this  study  are: 

1.  Flow  visualization  images  from  the  hovering  rotor 
capture  both  the  short-  and  long-wave  instabilities 
that  develop  in  the  near- wake  of  the  rotor  and  pre¬ 
cipitate  its  rapid  breakdown.  The  location  where 
the  wake  breaks  down  and  loses  its  structure  can 
be  seen  to  be  approximately  one  diameter  down¬ 
stream  of  the  rotor. 

2.  Flow  visualization  images  from  the  descending  ro¬ 
tor  operating  in  the  VRS  regime  show  the  merger 
of  the  individual  tip  vortices  and  the  formation 
of  a  thick  vortex  ring.  This  ring  would  remain 
just  above  the  rotor  plane  for  a  period  of  20-50 
revolutions  (depending  on  the  advance  and  de¬ 
scent  speeds,  and  on  the  rotor  rotation  rate)  before 
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abruptly  detaching  and  converting  away  upstream 
while  another  ring  would  begin  to  form. 

3.  Correlations  between  flow  visualization  images 
and  instantaneous  thrust  measurements  indicate 
a  severe  reduction  in  rotor  thrust  when  the  vortex 
ring  is  “attached  “  to  the  rotor,  followed  by  a  full 
recovery  of  thrust  once  it  is  shed. 

4.  The  regularity  of  the  vortex  ring  shed¬ 
ding/formation  process  over  100-500  revolution 
periods  indicates  that  the  process  is  quite  stable 
and  is  likely  dictated  by  the  size  of  the  ring  and 
the  amount  of  vorticity  it  can  contain. 

5.  Thrust  fluctuations  observed  in  the  VRS  regime 
were  most  severe  for  descent  angles  of  a=20-50° 
and  for  descent  speeds  of  V/14=1.0-1.5.  In  this 
region,  the  peak-to-peak  amplitudes  of  the  thrust 
fluctuations  were  approximately  80-95%  of  the 
mean  thrust.  Thrust  oscillations  of  this  magnitude 
can  have  disastrous  effects  on  the  performance  and 
control  of  a  rotorcraft  and  thus  deserve  further 
study.  A  major  question  that  should  be  consid¬ 
ered  concerns  whether  a  rotor  has  the  ability  to 
overcome  the  inflow  variations  imposed  on  it  by 
the  wake  oscillations. 

6.  In  VRS  cases,  the  vortex  ring  shedding  period  ex¬ 
hibits  a  decaying  power  law  dependence  on  the  ad¬ 
vance  speed,  indicating  that  the  ring  is  less  stable 
and  more  susceptible  to  being  swept  away  from  the 
rotor  as  advance  speed  increases.  Additionally,  the 
shedding  period  appears  to  increase  slightly  with 
collective  angle. 

7.  To  gain  a  more  complete  understanding  of  the  flow 
physics  in  this  flight  regime  it  will  be  necessary  to 
perform  PIV  experiments  on  the  critical  descent 
configurations  identified  in  this  paper.  By  collect¬ 
ing  quantitative  information  about  the  flow  field 
in  the  VRS  regime  it  should  be  possible  to  bet¬ 
ter  understand  the  nature  of  the  highly  regular 
vortex  ring  shedding/formation  process  was  ob¬ 
served  here.  And  with  knowledge  of  the  dynamics 
of  VRS,  it  should  eventually  be  possible  to  miti¬ 
gate  or  avoid  the  severe  thrust  fluctuations  asso¬ 
ciated  with  it. 
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